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The baryon acoustic oscillation (BAO) method for constraining the expansion history is adversely affected 
by non-linear structure formation, which washes out the correlation function peak created at decoupling. To 
increase the constraining power of low z BAO experiments, it has been proposed that one use the observed 
distribution of galaxies to "reconstruct" the acoustic peak. Recently Padmanabhan, White & Cohn provided an 
analytic formalism for understanding how reconstruction works within the context of Lagrangian perturbation 
theory. We extend that formalism to include the case of biased tracers of the mass and, because the quantitative 
validity of LPT is questionable, we investigate reconstruction in N-body simulations. We find that LPT does 
a good job of explaining the trends seen in simulations for both the mass and for biased tracers and comment 
upon the implications this has for reconstruction. 

PACS numbers: 



I. INTRODUCTION 

It has been known for many years that the coupling of pho- 
tons and baryons in the early universe results in significant 
features in the matter power spectrum tE 0, Ht]. Prior to re- 
combination, photons and baryons are tightly coupled and are 
well approximated by a fluid. Perturbations during this epoch 
do not grow, but instead excite sound waves which get frozen 
at recombination and manifest themselves as an almost har- 
monic series of peaks in the power spectrum, P{k), or equiv- 
alently a narrow feature in the correlation function, ^(r) (see 
10, Ht] for a detailed description of the physics in modern cos- 
mologies and 1 6] for a comparison of Fourier and configura- 
tion space pictures). These so-called "baryon acoustic oscil- 
lations" (BAO) can be used as a standard ruler to measure the 
expansion rate of the Universe, making the method an integral 
part of cuiTent and next-generation dark energy experiments. 

While the early Universe physics is linear and well under- 
stood, the low redshift observations are complicated by the 
non-linear evolution of matter, galaxy bias and redshift space 
distortions. The non-linear evolution leads to a coupling of k- 
modes and damping of the oscillations on small scales [4] and 
a small shift in their positions H IE S The damping 
of the linear power spectrum (or equivalently the smoothing 
of the correlation function) reduces the contrast of the feature 
and the precision with which the size of ruler may be mea- 
sured. 

In 1^ it was pointed out that much of the modification to 
the power spectrum comes from bulk flows and super-cluster 
formation. Since these large-scale flows are reasonably well 
measured by the survey, their effects can, in principle, be cor- 
rected. In lITll] a method was introduced for removing the 
non-linear degradation of the acoustic signature, sharpening 
the feature in configuration space or restoring/correcting the 
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higher k oscillations in Fourier space; this method has been 
tested on simulations by a number of groups fill [I^ [Htl- 
However, this method is inherently non-linear and therefore 
difficult to understand analytically. A study of this problem 
for the matter density using Lagrangian perturbation theory 
lfl4ll explained how the method "reconstructed" the BAO fea- 
ture, but also pointed out that it did not reconstruct the linear 
density field. We extend these results here - (i) generalizing 
the analytic theory to biased tracers, including explicit expres- 
sions for the reconstructed power spectrum to second order in 
the linear power spectrum, and (ii) testing the validity of the 
analytic expressions with a suite of N-body simulations. 

We compare the analytic theory to a set of 1024'^ parti- 
cle simulations run in periodic, cubical boxes of side length 
2/i^^Gpc with a TieePM code ifTsll . The simulations were 
initialized at z = 100 using second order Lagrangian pertur- 
bation theory, and the phase space information for all of the 
particles was dumped at z = 0, 0.3, 0.7 and 1.0. Multiple 
realizations, with different initial density fields, were run for 
each cosmology to reduce sampling effects (more details can 
be found in ifiol In addition to the dark matter particle 
data, halo catalogs were produced for each output using the 
friends-of-friends method iflTll with a linking length of 0.168 
times the mean inter-particle spacing. We work with halos 
above 10^"^ H^^Mq, i.e. containing more than 20 particles. 

We investigate one of the cosmologies considered in ifiol 
HI: ACDM, with = 0.25, = 0.04, h = 0.72, 
n = 0.97 and erg = 0.8. This is close to the current "best 
fit" cosmology and will serve as a realistic model to explore. 
Within this cosmology the acoustic peak in the correlation 
function is at ^ 110/i~^Mpc, with an intrinsic width set by 
the diffusion (Silk) damping scale of ~ 10 h^^Mpc. 

II. RECONSTRUCTION I: MATTER 

We begin our investigation of reconstruction by consider- 
ing the most conceptually simple case: reconstruction of the 
acoustic peak in the matter 2-point function. We start by re- 
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FIG. 1: The mass correlation function for our ACDM model at 2; = 
before (solid) and after reconstruction using a smoothing of i? = 
5 h^^Mpc (dotted) and R= 10 ft^^Mpc (short dashed). Non-linear 
evolution has partially erased the peak in the initial conditions (long- 
dashed) by z — 0, but it is somewhat restored by reconstruction. 

viewing the reconstruction algorithm of ifTlll and its interpre- 
tation within Lagrangian perturbation theory 13, and then 
compare its predictions with simulations. The following sec- 
tion extends this analysis to biased tracers. 

A. Algorithm 

The algorithm devised by [11] is straightforward to apply 
to a simulation and consists of the following steps: 

• Smooth the density field to filter out high k modes, 
which are difficult to model. 

• Compute the negative Zel'dovich displace- 
ment, s, from the smoothed density field: 
s(k) = — i(k/fc^)(5(k)iS(fc), where S is the smoothing 
kernel (see below). 

• Shift the original particles by s and compute the "dis- 
placed" density field, Sd- 

• Shift an initially spatially uniform distribution of parti- 
cles by s to form the "shifted" density field, 6s- 

• The reconstructed density field is defined as Sr = Sd — 
6s with power spectrum P,- (/c) cx ( | (5^ | ) . 

Following ifTHl we use a Gaussian smoothing of scale R, 
specifically 

We take advantage of the periodicity of the simulations to per- 
form all of these steps using fast Fourier transforms. The den- 
sity fields are constructed from the particle positions using a 
CIC assignment 1. 1 8.1 . 



Fig. [T] shows an example of reconstruction, for the ACDM 
model. By z = non-linear evolution has partially washed 
out the peak in the matter coiTelation function (solid line). 
However applying reconstruction with i? = 5 or 10 /i^^Mpc 
restores much of the original signal. Fig. |2] shows reconstruc- 
tion at the level of the density fields, for a thin shce through 
a piece of one of our simulations centered on a halo of mass 
4 X 10^'* h~^MQ. We see that reconstruction has 'reversed' 
the formation of collapsed structures, and yields a field that 
is visually similar to the initial density field. Note the final 
field has sharper, more pronounced peaks than either the ini- 
tial or reconstructed density fields, though the reconstructed 
field still has more prominent peaks than the initial field. 

B. Lagrangian perturbation theory 

Reconstruction naturally lends itself to a description in term 
of Lagrangian perturbation theory, which we briefly review 
here. The Lagrangian description of structure formation lfT9[ 
I20I I2TI1 relates the current (or Eulerian) position of a mass 
element, x, to its initial (or Lagrangian) position, q, through a 
displacement vector field \I/(q), 

x = q+*(q). (2) 

The displacements can be related to overdensities by ll22ll 

^(k) = J (fq e-^'^i (e-''' *(i' - 1) . (3) 

Analogous to Eulerian perturbation theory, LPT expands the 
displacement in powers of the linear density field, 5l, 

= + + • • • , (4) 

with ^(") being n'^ order in 5l- First order in LPT is equiv- 
alent to the well-known Zel'dovich approximation. 

In the simulations the rms (ID) displacement goes from 
6.1 h'^Wpc at z = to 3.9/i^^Mpc at z = 1, in excellent 
agreement with the expectations of the Zel'dovich approxi- 
mation. In fact the Zel'dovich rms displacements match those 
measured in the simulations at the percent level, better than we 
would expect given the size of the second order corrections. 

Using Eq. Q the power spectrum is 

P(fc) = j d?q e-^'^'i ((Je-*'=-'^*'(i)) - l) , (5) 

where q = qi — q2, and A\I/ — \I'(qi) — '4'(q2). Expand- 
ing the exponential in powers of ^ and using Eq. (lAlb re- 
produces the results of "standard" perturbation theory. How- 
ever, following [9], if we use the cumulant expansion theorem 
to expand the exponential and expand the resulting powers 
of k • A'S' using the binomial theorem we have two types of 
terms: those where the are all evaluated at the same point 
(which we can take to be the origin) and the rest. Leaving the 
first set of terms exponentiated while expanding the second 
set of terms in powers of yp, we find 

F(fc)=e-'='^'/2{Fi(fc) + ---} (6) 
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FIG. 2: A thin slice through a simulation showing the initial (left), reconstructed (middle) and final (right) density fields all smoothed with a 
Gaussian of 10 h~^Mpc. Each slice is centered on the (final) position of a halo of mass 4 x 10^'' h~^MQ. Note the final field has sharper, 
more pronounced peaks than either the initial or reconstructed density fields, though the reconstructed field still has more prominent peaks 
than the initial field. 
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TABLE I: The rms displacements of the "shifted" and "displaced" 
particles at 2 = and z = 1 as a function of the smoothing scale 
R. First order LPT correctly predicts the observed displacements (at 
the 10% level), with the agreement improving as the smoothing scale 
increases. 



with i[){x) — sin(a::)/a: the spherical Bessel function of order 
zero. This comparison also has the advantage of more clearly 
emphasizing the acoustic feature, which can be easily seen as 
a single peak in ^(r) at r ~ 110 /i^^Mpc. For presentation 
purposes we have smoothed all of the correlation functions by 
3 /i^^Mpc before plotting them - this reduces high frequency 
noise in the N-body simulations but has a minimal impact on 
the shape of the curves since this smoothing adds in quadra- 
ture to the ^ 10 /i^^Mpc intrinsic width of the features. Ob- 
servationally one could achieve similar effects by using broad 
but overlapping r bins. 

A second interesting statistic is the cross-spectrum between 
the linearly evolved initial field, 5^ and the fully evolved final 
field, (5/, 



where PL{k) is the linear theory power spectrum, S is pro- 
portional to the rms Zel'dovich displacement (i.e. final minus 
initial particle positions to linear order) 
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and explicit expressions for the higher order terms may be 
found in Jgt] and Appendix |A] The exponential prefactor de- 
scribes the broadening of the acoustic peak seen in Fig. [T] 
some of the additional terms lead to a slight change in the 
peak position |[l,[l3>li3l- The rms displacement of an individ- 
ual particle is E/\/2. 

The effects of the exponential prefactor are most easily seen 
by considering the correlation function. Furthermore, La- 
grangian perturbation theory, like several other perturbation 
theory schemes, performs better at predicting the large-scale 
coiTelation function than the power spectrum, since it fails 
to accurate ly p redict broad-band power which contributes at 
small r ^ uM- For these reasons, we shall present most of 
our comparisons between theory and simulation in configura- 
tion space, i.e. we shall present 

P{k)j,{kr) = j - A^{k)Mkr) (8) 



Gf{k) 



SL{k)5)(k) 



(9) 



sometimes referred to as the propagator ||8|]. The relevant 
physics in this case is more cleanly visualized in Fourier 
space, since it shows the decorrelation between the initial field 
and the processed field which becomes a convolution in con- 
figuration space. Fits to numerical simulations [6] and a va- 
riety of analytic arguments [6, 7, 8, 9, JJD, J4, 23J, including 
Lagrangian perturbation theory, suggest that 



Gf{k) 



-(fcE)V4 



(10) 



i.e. that the damping is half as strong as in the power spectrum 
(see AppendixlAlfor expressions beyond leading order). 

It is straightforward to repeat these steps for the recon- 
structed field [14J. We assume that the density field is 
smoothed on a large enough scale that s can be approximated 
as s = — i(k/fc^)(5L(k)5(fc). We can then compare the three 
contributions to the power spectrum {Pss, Pdd and Psd) to get 
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the reconstructed power spectrum ifTill 

+2e-'='^'^/25(fc)5(fc) (11) 
+ e-'='^'^/252(fc)}Pi(fc) + ... 

where S = 1 — 5, and as before, the higher order terms are 
AppendixlAl There are now three smoothing terms, {Y,ss, ^dd 
and J^sd) defined by 

^1 = ^ J dqPL{p)S'{p), (12) 
^Id^^ I dqPL{p)S'{p), (13) 

and = (S^ + S^d)/2 (see Table ID- As pointed out in 
Ref. [14^, all of these smoothing scales are smaller than the 
nonlinear smoothing E, explaining why the acoustic feature 
is sharpened after reconstruction. A related calculation (see 
Appendix [All yields the propagators 



Gf = e-''^^^^^ + ... (14) 

Gd = e-^'^'^^/^S + ... (15) 

- e-'^'^-/^ [-S] + ... (16) 

Gr = Gd — Gs , (17) 



with the higher order terms in the Appendix. 

C. Comparison with simulations 

To begin we compare the predictions of perturbation the- 
ory for the propagator to calculations of the same quantity in 
N-body simulations. This isolates the damping behavior from 
the mode-coupling |10]. Fig.|3]shows the different contribu- 
tions to the reconstructed propagator The theoretical predic- 
tions for Gf are in reasonably good agreement with the re- 
sults, with the theory showing slightly weaker damping than 
the simulations. (Small changes to the theoretically predicted 
E can bring the results into much better agreement, but we 
will not make such ad hoc changes here.) The agreement is 
somewhat worse for some pieces of the reconstructed propa- 
gator. In particular the simulations show that the reconstructed 
field retains better memory of its initial conditions (Gr ~ 1) at 
intermediate scales than LPT predicts, with perturbation the- 
ory giving too much power at high k. The over-prediction at 
high k is not of particular concern, since at these scales the di- 
mensionless power exceeds unity and we would expect pertur- 
bation theory to be breaking down. Out to fc ~ 0.2 Mpc^^, 
where ~ 1, perturbation theory agrees with the simula- 
tions at the better than 10% level! We emphasize that this 
level of agreement comes from the inclusion of the 2""^ order 
contributions, with the dominant correction coming from the 
i?i term (see Appendix [All. 

Fig. |4] shows the corresponding figure for the correlation 
functions, broken down into the same components. Note the 
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FIG. 3: The cross-correlation between the linearly evolved initial 
field and the fully evolved final field, displaced field, shifted field 
and the reconstructed field for ACDM at z = (see text). The 
points show the results of N-body simulations while the lines show 
the predictions from Lagrangian perturbation theory 1 14]. The solid 
line and diamonds represent G/(fe), the dotted line and circles repre- 
sent Gd(k), long-dashed line and crosses represent Gs{k) and short- 
dashed line and triangles represent Grik). We have used a smooth- 
ing ofR= 10/i"^Mpc. 

excellent agreement for the displaced and shifted fields, but 
less good agreement for the final and reconstructed fields. In 
this figure the level of agreement between and the theory is 
worse than the comparable figure in fl^. This is most likely 
due to the lower redshift and different cosmology we have 
chosen (see also 1 16]). The sense of the disagreement in both 
^/ and is the same however, indicating that the Lagrangian 
perturbation theory of reconstruction is working better in a 
differential than absolute sense. As above, a small change in 
the relevant E could slightly improve the agreement with sim- 
ulations, which may argue for leaving E as a free parameter 
when fitting to data. We will not pursue such modifications 
further here. 



III. RECONSTRUCTION II: BIASED TRACERS 

Unfortunately we don't directly measure the mass field in 
galaxy surveys, we measure the distribution of biased tracers. 
Here we investigate how the biasing of the tracers affects re- 
construction. Rather than attempt a 'realistic' galaxy model, 
we shall concentrate on mass limited samples of halos when 
comparing LPT to the simulations. None of the essential as- 
pects are lost with this simplification. 

Reconstruction assumes that we can estimate the appro- 
priate shifts from our smoothed, biased, density field. This 
requires that the smoothed halo field be a multiple of the 
smoothed mass field with known constant of proportionality 
(the bias). In the simulations we estimate the bias from the 
k ~ Q limit of the propagator, in observations it would need 
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FIG. 4: The correlation functions of the fully evolved final, dis- 
placed, shifted, and reconstructed mass fields for ACDM at z — 
(see text). As in Fig. [3] the points show the results of N-body sim- 
ulations while the lines show the predictions from Lagrangian per- 
turbation theory [140. Solid line and diamonds represent dotted 
line and circles represent ^d, long-dashed line and crosses represent 
^s, and short-dashed line and triangles represent ^r- We have used a 
smoothing of R — 10 h^^ Mpc. 



to be determined in a different manner. 

If we keep the denominator in Eq. (|9|l as the linear mass 
power spectrum, the lowest order modification to the prop- 
agator is to multiply by the linear bias of the tracer The 
gross shape of G{k) is unaltered, since the exponential damp- 
ing is unchanged, being generated by the velocities which are 
sourced by the mass field not the halo field. At higher order, 
the cross terms between the linear and rtth order terms are 
modified and introduce an additional dependence on the bias 
||24I1 . as shown in Fig.|5] LPT predicts that the halo propagator 
falls slightly more slowly to high k than the mass propagator 
and the decline is slower the higher the mass threshold. This 
means that the halo propagator departs more from the Gaus- 
sian form than the mass propagator It is possible that this is 
related to the special locations in the velocity field that rare, 
highly biased peaks occupy (e.g. |'25','26'l). However, the dif- 
ference is small, as shown explicitly in Fig.|5] 

Fig. |6] shows the different propagators for halos more mas- 
sive than 10^^ H'^Mq in the simulations and in theory. As 
was the case for the mass, the asymptote at high k is not well 
determined by the theory but the agreement at low k is quite 
good. Gs is the same as for the mass, and again the agreement 
between simulation and theory is good. The match between 
simulations and theory for G/ is quite good. Perturbation the- 
ory is correctly predicting the low k asymptote of Gd, which 
is no longer zero but 6—1, though it doesn't match the shape 
as well as for the mass. Once more the N-body simulations 
predict a which increases slightly at intermediate k and is 
above the theory for fc ~ 0.1 - 0.2 /iMpc"^ 

To lowest order (see Appendix |A] for 2"^^ order contribu- 
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FIG. 5: The cross-correlation between the linearly evolved initial 
field and the fully evolved final field for the mass (dashed line and 
circles) and for halos above 10^^ H^^Mq (solid line and squares) 
in ACDM at z = 0. The dotted line shows the mass propagator 
multiplied by b ~ 1.6. 
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FIG. 6: The cross-correlation between the linearly evolved initial 
field and the evolved, displaced, shifted and reconstructed fields for 
halos above 10^"^ h~^MQ. Diamonds and the solid line show the 
final field, crosses and the long-dashed line the shifted field, circles 
and the dotted line the displaced field and the triangles and short- 
dashed line the reconstructed field. 

tions) the reconstructed field has 

Pi"^(fc) =Pi(fc){e-'='^==/252(fc) 

+2e-'='^'^/2 [S{k)S{k) + {b- l)5(fc)] 
+ e-k'^lJ^ [S^(^k) + 2(6 - l)S{k) + (6 - 1)2] I 

(18) 

which reduces to Eq. ( fTTl i in the limit 6 ^ 1. Note that 

P!:°\k) ^ 62Pi(fc) as fc 0, as expected, and Pi"'(fc) -> 
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FIG. 7: The correlation functions for tiie evolved, displaced, shifted 
and reconstructed fields for halos above 10^'^ h^^AlQ. Diamonds 
and the solid line show the final field, crosses and the long-dashed 
line the shifted field, circles and the dotted line the displaced field 
and the triangles and short-dashed line the reconstructed field. 

62Pi(fc) exp[-fc2i;2/2] in the limit that = Y.dd = ^sd- 

Fig. |2] shows how well this expression, plus the 2"^^ order 
contributions, matches the simulations. As with Fig. H] the 
agreement is overall quite good, slightly better than for the 
mass in the case of ^/ and i^^- As in that case, a slight increase 
in the S can improve the agreement somewhat, but we have 
left the theoretical predictions unchanged. 

Just as with the matter field, the smearing of the acoustic 
peak is reduced by reconstruction. In fact there is relatively 
little difference between the biased and unbiased tracers in 
this respect. 

IV. CHANGE IN THE PEAK LOCATION 

The above sections demonstrate that the LPT provides a 
good description of how reconstruction reduces the smooth- 
ing of the acoustic feature, both for the dark matter and ha- 
los. Recent simulations |7] have also found that reconstruc- 
tion corrects the ~ 0.5% change in the acoustic scale caused 
by nonlinear evolution. It is therefore interesting to see how 
this is manifest within Lagrangian perturbation theory. 

In perturbation theory the change in the acoustic peak loca- 
tion comes about because there are second-order corrections 
to P{k) which are out-of-phase with the linear theory oscil- 
lations 1 6, 7, 8, 9, 14]. The out-of-phase component is quite 
similar to the derivative of Pl so, by Taylor's theorem, this 
addition is akin to a change in the characteristic frequency of 
the oscillation. We consider the analogous terms for the re- 
constructed power spectrum below, in order to explain how 
reconstruction suppresses such changes. 

These out-of-phase components come about because of the 
structure of the mode-coupling terms (the Qn in the notation 
of Appendix lAl and Ref. |24]), and this structure is modified 
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FIG. 8: The out-of-phase pieces of the power spectrum of halos more 
massive than 10^^ H^^Mq as predicted by perturbation theory. To 
emphasize the oscillations, each spectrum has been divided by the 
"no wiggle" form of Ref. |5] and has had a 4"^ order polynomial (in 
k) subtracted. The dotted line shows the linear theory (divided by 
2). The solid line is the out-of-phase or mode-coupling pieces of Pf, 
which can be compared to (IPl / din k (long-dashed line) 1 10]. The 
short-dashed line shows that reconstruction reduces the amplitude of 
the out-of-phase terms and hence the change in the location of the 
acoustic peak in ^(r). 

by reconstruction in such a way as to reduce the amplitude of 
the out-of-phase contribution |14]. Figure [8] shows the out- 
of-phase terms, with the broad-band shape removed to focus 
on the oscillatory structure, compared to the in-phase acous- 
tic signature in the linear theory. Note that the modification 
of the mode-coupling terms detailed in the Appendix drasti- 
cally reduces the amplitude of the out-of-phase terms in the 
reconstructed spectrum, and hence the change in the acoustic 
scale. This explains why the change in the peak location seen 
in simulations is reduced by reconstruction. 

V. DISCUSSION 

Acoustic oscillations in the photon-baryon fluid prior to de- 
coupling leave an imprint both in the cosmic microwave back- 
ground anisotropy power spectrum and the matter power spec- 
trum. A comparison of these features at different redshifts 
provides one of the most promising routes to constraining 
the expansion history of the Universe. Unfortunately at low 
redshift, where the accelerated expansion of the Universe is 
strongest, non-linearities wash out much of the acoustic infor- 
mation. 

Recently the authors of Ref. ifTlll proposed a method for re- 
covering much of the lost information, or reconstructing the 
acoustic peak. Unfortunately, the method is inherently non- 
linear and therefore difficult to understand analytically. A 
study of this problem in Lagrangian perturbation theory 13, 
for the mass field, shed some light on how the algorithm re- 
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suited in tighter constraints on the acoustic scale, but the quan- 
titative validit y of Lagrangian perturbation theory is question- 
able (see e.g. udH for a recent survey) and we typically study 
biased tracers of the mass. 

We have validated and extended the analytic insights devel- 
oped in 131, computing a variety of statistics of both the mass 
density field and the dark matter halo density field using La- 
grangian perturbation theory which we then we then compare 
to the same quantities measured in a large suite of N-body 
simulations. 

As emphasized in ifTill . reconstruction does not generate 
the initial power spectrum or correlation function, but it does 
serve to sharpen the peak and reduce the change in the peak 
location associated with non-linearity. We demonstrate ex- 
plicitly that both of these points remain true for biased trac- 
ers. The amount by which the non-linear smearing is reduced 
is comparable for biased tracers and for the mass, since it is 
generated by bulk flows which are sourced by the mass den- 
sity independent of the form of the tracer. The fact that peaks 
form in special locations in the density field appears to have a 
very small effect. The reduction in the peak location change 
due to reconstruction is at least as dramatic for biased tracers 
as for the mass, with the out-of-phase component responsible 
for the change being reduced in amplitude by the process of 
reconstruction. A discussion of by how much the peak posi- 
tion changes depends on a detailed description of the fitting 
methodology and the sample under consideration, but if we 
model the observed spectra as in ifToll we find that reconstruc- 
tion reduces the position change by a factor of 2 — 4 for mod- 
erately biased tracers like those investigated here. 

We conclude that Lagrangian perturbation theory, while not 
perfect, provides a good framework for thinking about recon- 
struction. It explains in a natural way how reconstruction 
works, and how it achieves a reduction in the smearing and 
position of the acoustic peak generated by non-linear evolu- 
tion. The predictions of LPT agree to within several percent 
with the results of N-body simulations on the large scales most 
relevant to acoustic oscillations, for both biased and unbiased 
tracers. While not shown explicitly in figures here, perturba- 
tion theory becomes an increasingly good description of the 
simulations at higher redshift, though the need for reconstruc- 
tion beyond z ~ 1 is greatly reduced. 
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APPENDIX A: BEYOND LEADING ORDER 

Lagrangian perturbation theory allows us to compute cor- 
rections to the lowest order expressions for P{k) and G{k) 
listed in the text. Here we give the 2"^^ order contributions, 
following ll9l l24ll and ifToll . The notation and procedure is bor- 
rowed heavily from these works, to which we refer the reader 
for more details. 

Recall that the displacement is expanded in powers of the 
linear density contrast, (5^, as iIztIi 



*(")(k) 



n 



(27r)3 



X L(") (ki , • • • , k„ , k)<5i (ki ) • • • Sl (k„ ) (Al) 

where the L'"^ have closed form expressions, generated by 
recurrence relations. For example. 



(A2) 



is the well known Zel'dovich displacement, which is 1'^* order 
LPT. 

The density field for a biased tracer can be defined by the 
displacement field ^ (q) and a function of the smoothed initial 
density field in Lagrangian space, F[Sl{(1)], 

Jobj(x) = J d\F[SL{q)]S^^\^ - q - *) , (A3) 

where x and q are the Eulerian and Lagrangian positions and 

(3) 

6)2, is the 3D Dirac 6 function. The power spectrum for such 
tracers can then be written as [9, 24] 

dXi d\2 



P{k) = j d^qe-"^"^ j 



27r 27r 

,i(Ai5i,(qi)+A2<5t(q2))+ik[*(qi)-*(q2)] 



F{\i)F{X2)x 



, (A4) 



where q = qi — q2 and F is the Fourier transform of F. 
The distribution-averaged derivatives of F{X), {F') and {F"), 
characterize the bias of the sample under consideration. Ex- 
pressions for the case of peaks in the initial density field 
(i.e. peaks bias) can be found in ll24ll . For the halos consid- 
ered in flie text (M > 10" H'^Mq) we have {F') = 0.55 
and {F") = -0.37, with large-scale bias 1.55. 

To obtain the propagator we cross-correlate Eq. ( IA3I ) with 
a field defined by cxp{iXSL)', is then simply obtained by 
taking the rt-th derivative with respect to A and setting A to 
zero iflOll . This allows us to follow a procedure similar to that 
in Eq. ( IA4I ). 

The algebra now follows through as in 1^ |23l using the 
cumulant expansion theorem, and coflecting all zero-lag cor- 
relators to yield, e.g. 



{SlSo 



hi) 



qe 



-ikq 



(A5) 
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where we have omitted the exponential damping terms for 
brevity and defined l.24il 



mim2 — V J 



X ([<5i(qi)]"i[5L(qi)]"^[k*(qi)]™i[k*(q2)]™^)c , 

(A6) 

with (• • ■)c denoting the connected moments. 
Straightforward algebra then yields 



where S = (1—5) and is evaluated using P^S inside 
the integral Eq. ( IA8I ). 

The power spectrum can be evaluated in a similar fashion, 
with the three contributions being 



and 
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where \2i 



and 
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+ {F')[Pl + -{Ri+R2} 
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dr PL{kr)Rn{r) (A8) 
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(A16) 
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while for the power spectrum, omitting the damping terms, 

Pobj (X (l + (0)'^L + ^Ql + yQ2 + iQ3 



(F" 



+ (F')'[Q9 + Qn]+2(F')(F")gi2 
+ i(n'Qi3 + y(l + (n)'[i?i+i?2] 



with 



21 



Qn{k) 



{^ + {F'))Ri 



(All) 



(2^)^ .0 
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dr PL{kr) j d^j. 



Pl 



Q„(r,Ai) (A12) 



and expressions for the Qn can be found in f2A 

Extending these results to reconstruction is now relatively 
straightforward, assuming that the smoothed density field can 
be well approximated by the linear field. The shifted field 
has = 'i/'^'iS and no higher order contributions, while the 
displaced field can be obtained from 'S' with the replacement 
^(1) _> ijr(i) [1 _ s] with *(">2) unchanged. This yields 



(SlS,) cx -PlS 
ihSd) « P^S+^R, 



-R^ 



(d) 



(F') [Pl + -{Ri+ R2} ] (A14) 
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2(F')(F")Ql2"''' 



(A17) 



where we have again omitted the damping terms and the su- 
perscripts indicate which Pl are to be replaced with PlS, 
PlS etc. For the Qn there are 4 possible smoothing terms, 
and we have indicated no smoothing with a 1, 5 with an s and 
S with a d. The first two terms have argument kr and the sec- 



ond have argument k ^/l 
example 



^{Isds) 
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2r^ in Eq. dAHb . Thus for 
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(A 13) with y = ^Jl + r"^ ~ 2r/i. 
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